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ABSTRACT 


A  PRN  (packet  radio  network)  is  a  collection  of  geographically 
distributed,  possibly  mobile  users  where  each  user  is  capable  of  trans¬ 
mitting  and  receiving  messages  over  a  shared  broadcast  medium.  In  a 
PRN,  messages  are  divided  into  packets,  which  may  be  fixed  or  variable 
in  length,  and  each  packet  is  transmitted  through  the  network  individually. 
Packets  are  assembled  at  their  destinations  to  reconstruct  the  original 
messages. 

.  The  data  traffic  in  a  PRN  is  characterized  by  specifying  the  average 
message  arrival  rates  to  the  network  for  each  o-d  (origin-destination) 
pair.  A  set  of  o-d  rates  is  called  feasible  if  th're  exist  network 
protocols  under  which  the  number  of  packets  in  the  network  still  not 
delivered  to  their  destinations  remains  finite  with  probability  one. 

The  capacity  region  of  a  PRN  is  defined  to  be  the  set  of  all  feasible 
sets  of  o-d  rates. 

In  this  thesis,  PRNS:  are  studied  from  the  viewpoint  of  feasibility, 
i.e. ,  we  take  an  arbitrary  set  of  message  input  rates  as  given  and  try 
to  determine  if  it  is  feasible.  Our  main  conclusion  is  that,  unless 
P  “  NP»  there  exists  no  practical  algorithm  for  characterizing  the 
capacity  region  of  a  PRN,  in  the  sense  that  the  decision  problem  regarding 
the  feasibility  of  a  given  set  of  o-d  rates  is  NP-conplete. 
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CHAPTER  I.  INTRODUCTION 

1.1  Background  Information 

Recently  there  have  been  attempts  to  extend  the  domain  of  data 
communications  to  networks  of  geographically  distributed  mobile  radio 
users.  Radio  networks  have  traditionally  been  used  for  voice  communi¬ 
cations  such  as  in  police  cars,  emergency  vehicles,  etc..  We  isider 
here  radio  networks  designed  mainly  for  computer  data  communice  t»®. 

Such  networks  present  a  multitude  of  problems  which  must  be  re  to 

make  the  idea  a  feasible  one. 

One  of  the  problems  results  from  the  bursty  nature  of  computer 
data  traffic.  Typically,  a  long  period  of  silence  is  followed  by  a 
sudden  burst  of  huge  amounts  of  data  which  must  be  transmitted  through  the 
network  to  the  desired  destination. 

The  fact  that  a  station  need  not  be  (and  typically  is  not)  within 
the  transmission  range  of  every  other  station  further  complicates  the 
matter.  Unlike  a  wire  network,  transmissions  may  interfere  with  one 
another  and  this  causes  failure  to  detect  messages  correctly. 

It  is  a  challenging  task  to  design  network  protocols  which  are 
distributed  in  nature  and  which  satisfy  the  service  demands  of  the 
users.  In  this  thesis  we  study  two  well-known  protocols,  namely 
ALOHA  and  TDMA,  under  a  model  which  is  refined  of  all  non-essential 
features  so  as  to  simplify  the  analysis. 

For  a  general  survey  of  packet  radio  networks  we  refer  the  reader 
to  [Hi  and  [KiSBK  78] . 


1.2  The  PRN  Model 


We  assume  that  messages  to  be  transmitted  from  one  station  to 
another  are  first  framed  into  packets  which  may  be  of  fixed  or  variable 
length.  Subsequently,  each  packet  is  transmitted  separately.  The 
packets  contain  the  information  necessary  for  their  travel  through  the 
network,  such  as  the  identities  of  the  origin  node,  the  destination  node 
etc.  . 

We  shall  represent  a  PRN  (packet  radio  network)  by  a  directed  graph 
G  »  (N,A).  (For  definitions  of  graph  theoretic  terms  used  in  this 
thesis,  see  e.g.  [BJ.)  To  each  station  in  the  network  there  corresponds 
a  node  of  G,  and  conversely.  Throughout  the  thesis  we  use  the  words 
user,  node,  and  station  interchangeably. 

For  any  two  distinct  nodes  a,  b  e  N,  there  is  a  link  (a,b)  iff 
node  b  is  within  the  transmission  range  of  node  a.  We  do  not  assume 
that  (a,b)  6  A  implies  (b,a)  €  A. 

We  shall  study  only  PRN's  with  finite  numbers  of  users;  accordingly 
G  -  (N,A)  will  always  be  a  finite  graph. 

The  graphical  representation  of  a  PRN  will  be  simple,  that  is, 

G  ■  (N,A)  will  not  contain  self-loops  or  multiple  links  from  one  node 
to  another. 

Throughout  the  thesis  L  will  denote  the  number  of  links  in 
G  -  (N,A) ,  i.e.  L  =  |a|. 

A  PRN  is  said  to  be  connected  if  packets  can  be  routed  from  every 
user  to  every  other  user.  Equivalently,  a  PRN  is  connected  if  the 
corresponding  graph  G  is  connected,  i.e.  if  there  exists  a  directed 
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path  from  every  node  to  every  other  node  in  G. 

In  our  model/  each  user  of  the  PRN  has  a  single  transmitter  and  a 
single  receiver,  thus,  a  user  can  transmit  or  receive  only  one  packet 
at  a  time. 

When  a  user  transmits  a  packet,  it  is  possible  that  all  stations 
within  its  transmission  range  receive  this  packet  without  error;  however, 
the  packet  may  be  intended  only  for  a  certain  subset  of  them.  In  fact, 
in  this  thesis  we  shall  assume  that  each  transmitted  packet  has  exactly 
one  intended  receiver. 

If  a  packet  is  received  without  error  by  its  intended  receiver, 
the  transmission  is  said  to  be  successful.  We  assume  that  the  only 
source  of  unsuccessful  packet  transmissions  is  interference.  A  packet 
which  has  not  been  received  successfully  is  said  to  have  suffered  a 
collision. 

Let  (a,b)  be  a  link;  we  say  that  a  packet  is  transmitted  over  link 
(a,b)  if  the  packet  is  transmitted  by  node  a  with  the  intended 
receiver  being  node  b.  We  say  that  link  (c,d)  conflicts  with  link 
(a,b)  if  simultaneous  transmissions  on  links  (c,d)  and  (a,b)  either  are 
precluded  (when  c**a,  dj*b)  or  cause  a  collision  at  node  b  (when  c*b,  or 

(c,b)  is  a  link) .  We  define  to  be  the  set  of  links  which  conflict 
with  (a,b) .  Therefore,  we  have  that 

-  {(c,d)  e  A:  c  ^  a,  (c,b)  8  A)  U 
{ (a,d)  8  A:  d  *  b}  U  {(b,d)  e  A> 

Sometimes  we  shall  consider  the  set 

'  CIb  "  cab  “  ^b'd>  S  A)  U  { (c,d)  e  A  t  (c,b)  8  A> 
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Remarks: 

1)  <c,d)  SCil  ^  (.,b)  e  Ccd 

*  * 

2)  for  any  two  links  (a,b),  (c,b)  e  A 

This  completes  the  description  of  the  physical  structure  of  the 
network.  In  the  next  section,  we  shall  look  at  another  aspect  of  the 
network;  namely,  the  data  traffic. 


1.3  Data  Traffic  in  a  PRN 

The  origin  of  a  packet  is  the  node  at  which  the  packet  enters  the 
network;  the  destination  of  a  packet  is  the  node  at  which  the  packet 
leaves  the  network.  In  our  model  each  packet  has  a  unique  destination. 

With  every  o-d  (origin-destination)  pair  we  associate  two  random 
processes  parametrized  on  the  time  interval  [0,°° ).  The  x-y  arrival 
process  (t)  is  the  number  of  packets  which  have  arrived  at  origin  x 
with  destination  y  in  the  interval  [0,tj,  t  >  0.  The  x-y  departure 
process  (t)  is  the  number  of  packets  which  have  been  delivered  from 
origin  x  to  destination  y  in  the  time  interval  f0,t),t^  0.  We 
assume  that  prior  to  time  0,  there  are  no  packets  in  the  network. 

Let  W  be  the  set  of  all  o-d  pairs;  let  |w|  “  w ;  and  let  the  o-d 
pairs  be  labeled  by  integers  l,...,w. 

we  collect  the  arrival  and  departure  processes  associated  with 
all  o-d  pairs  into  w-dimensional  vectors  R(t)  and  D(t),  respectively. 
Notice  that,  R(t)  ^  D(t)  and  1  •  (R(t)  -  D(t))  is  the  number  of 
packets  in  the  network  at  time  t. 


The  process  D(t)  depends  both  on  R(t)  and  the  decisions  made  by 
all  the  protocols  that  determine  the  flow  of  packets  through  the  net¬ 
work. 

The  primary  objective  in  a  PRN  is  to  find  network  protocols  such 

that 

lim  lim  p{T  •  (R(t)  -  D  (t) )  <  a}  ■  1  •  (1.1) 

a** 

In  this  thesis,  we  restrict  our  attention  to  the  following  class 
of  o-d  arrival  processes. 

1)  R(t)  is  independent  of  the  network  protocols. 

2)  For  all  iew,  there  exists  r^  :>  0  j  and,  given  - 

1  • 

e  >  0  ,  $  >  0  we  can  find  T(e,<S)  >  0  such  that 

iR.(t)  . 

p{|— t  -  rj  >  e  }  <  «  vt  >  T(e,6) . 

Clearly ,  r^  with  the  above  property  is  unique  when  it  exists.  We 
define  r  to  be  the  column  w-vector  whose  i^1  row  equals  r^j  r  is 
called  the  mean  o-d  (arrival)  rate  vector. 

The  mean  o-d  rate  vector  r  will  be  used  to  characterize  the 
arrival  process  R(t) .  We  shall  say  that  r  is  feasible  if  there 
exists  network  protocols  such  that  (1.1)  is  satisfied. 

We  define  the  capacity  region  of  G  ■  (N,A)  to  be  the  following  set: 

C(G)  *  {r  i  r  is  feasible) 


It  is  of  extreme  interest  to  determine  the. region  C(G).  A  second 
problem,  which  is  also  interesting,  is  to  determine  whether  a  particular 
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mean  o-d  rate  vector  is  feasible  under  a  particular  set  of  network 
protocols.  By  network  protocols  we  essentially  mean  the  routing  and 
link- level  multi-access  schemes.  The  interaction  between  routing  and 
access  schemes  makes  it  difficult  to  study  the  feasibility  problems. 
Therefore^  we  examine  multi-access  schemes  in  isolation  from  routing 
schemes  by  assuming  that: 

1)  r  -  0  if  (x,y)0  A 

2)  Every  (x-y)  packet,  i.e.  a  packet  with  origin  x  and 
destination  y,  is  transmitted  directly  over  link  (x,y) 

When  we  consider  the  feasibility  problem  under  these  assumptions 

we  shall  use  the  column  L-vector  f  instead  of  r.  The  mean  link 

(arrival)  rate  vector  f  has  one  element  for  each  link  (a,b)  e  A, 

which  equals  r^  .  Thus,  f  is  the  mean  link  rate  of  (a,b) .  Under 

the  above  assumptions ,  the  analogues  of  R  (t) ,  R(t) ,  D  (t)  and 

xy  xy 

D(t)  will  be  denoted  by  F  (t) ,  F(t),  S  (t)  and  S(t),  respectively. 

xy  xy 

For  example,  F  (t)  -  S  (t)  is  the  number  of  packets  waiting  at 
xy  xy 

node  x  to  be  delivered  to  node  y  over  the  link  (x,y)  at  t 
(including  a  packet  in  transmission ,  if  any). 

With  this  notation,  we  can  now  define  the  two  main  problems 
we  shall  study  in  this  thesis. 

Definition . 

The  r  -  feasibility  problem 

Given  a  PRN  G  -  (M,A)  and  a  w-vector  r  >  is  r  feasible,  i.e., 
does  r  belong  to  C(G)7 
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Definition. 

The  £  -  feasibility  problem 

Given  a  PRN  G  *  (N,A) ,  a  multi-access  scheme ,  and  an  L-vector  f  _>  0, 
is  ?  feasible  under  the  given  multi-access  scheme? 

We  shall  study  the  f-feasibility  problem  under  two  extreme  ways  of 
multi-accessing/  namely  ALOHA  and  TDHA.  The  r-feasibility  problem  will 
be  studied  in  connection  with  TDMA.  Before  we  give  an  outline  of  our 
results  on  the  feasibility  problems ,  we  shall  review  some  of  the  relevant 
work  in  this  area. 

1.4  Review  of  Earlier  Work 

In  a'B.S.  thesis  Shapiro  (S791  analyzes  an  Aloha  scheme  for  PRNs 
and  obtains  conditions  for  minimizing  average  system  delay  in  terms  of 
routing  variables. 

A  different  approach  to  the  analysis  of  Aloha  type  systems  is 
taken  by  Sidi  and  Segall  [SS(1)81]  who  consider  the  PRN  as  a  network 
of  interfering  queues  and  obtain  approximate  results  for  the  steady- 
state  packet  distribution  in  the  buffers  of  the  PRN. 

PRN  a  have  also  been  studied  under  different  assumptions  and 
multi-access  schemes,  e.g.  [SS(2)81]. 

1.5  Thesis  Outline 

In  Chapter  2,  the  1-feasibility  problem  is  studied  under  slotted 
Aloha.  We  examine  conditions  <-or  stability  from  a  feasibility  view¬ 
point;  routing  variables  and  link  delays  do  not  enter  into  the  analysis. 
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In  Chapter  3,  we  formulate  and  state  a  number  of  facts  about 
TDMA  schemes,  among  these  facts  are  the  NP-completeness  of  the 
f-feasibility  problem  under  TDMA,and  the  NP-completeness  of  the 
r-feasibility  problem. 

The  fourth  chapter  consists  of  NP-completeness  proofs  of  the 
feasibility  problems  and  an  alternative  formulation  of  TDMA  schemes. 

In  Chapter  5,  we  show  why  not  even  a  polynomial-time  approximation 
algorithm  is  likely  to  be  found  for  the  NP-complete  feasibility 
problems. 

The  main  result  of  this  thesis  is  the  conclusion  that,  unless 
P  »  NP,  there  exists  no  practical  algorithm  for  characterizing  the 
capacity  region  of  a  PHN. 
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CHAPTgR  II.  SLOTTED  ALOHA  1  ti  PENS 
2.1  Introduction 

In  this  chapter  we  study  the  f-feasibility  problem  under  slotted 

i 

Aloha.  We  assume  familiarity  with  Aloha  as  it  is  used  in  single¬ 
receiver  multi-access  communication  systems.  The  scheme  considered  here 
is  a  simple  extension  of  the  well-known  Aloha  systems;  we  do  not  claim 
that  it  is  original  or  practical.  Our  purpose  is  to  study  this  scheme 
along  the  lines  outlined  in  Chapter  I. 

The  main  difficulty  in  the  study  of  Aloha  in  PENS  is  encountered 
at  the  modelling  stage.  Our  model  must  account  for  the  essential 
features  of  Aloha  and  yet  it  must  be  tractable.  We  can  make  the  model 
more  and  more  sophisticated  only  to  render  it  useless.  On  the  other 
hand,  an  oversimplified  model  may  lead  to  erroneous  conclusions.  The 
modelling  issue  is  further  complicated  by  the  variability  of  network 
topology  and  data  traffic. 

Our  analysis  of  Aloha  is  based  on  a  model  which  is  as  simple  as 
it  can  be.  We  introduce  the  equilibrium  hypothesis  and  discuss  the 
meaning  of  equilibrium.  The  analysis  provides  results  which  on  the 
whole  conform  with  intuition  provided  by  the  single-receiver  Aloha 
systems.  However,  some  of  the  results  are  non -obvious  and  seem  to  be 


useful. 


2.2  The  Aloha  Model 


We  consider  a  slotted  Aloha  scheme  where  all  stations  are  synchronized 
so  that  packet  transmissions  occur  only  in  globally  defined  time  slots. 
Packets  are  fixed  in  length  and  the  duration  of  a  slot  is  long  enough  to 
accommodate  the  time  it  takes  to  transmit  a  packet  plus  any  delays 
associated  with  propagation  and  detection  of  packets .  We  assume  pro¬ 
pagation  delay  is  negligible  relative  to  the  length  of  a  slot,  but  we  do 
not  ignore  it  altogether. 

In  Aloha,  transmitted  packets  may  suffer  collisions,  and  when  they 
do,  they  must  be  retransmitted.  Thus,  there  are  three  random  processes 
associated  with  each  link  that  we  have  to  distinguish: 

1)  The  process  of  new  packet  arrivals, 

2)  The  process  of  successful  packet  transmissions, 

3)  The  scheduling  process  of  packet  transmissions. 


Following  the  notation  introduced  in  Chapter  1  in  connection  with 


the  f-feasibility  problem,  we  define  F^tn)  n  6  {l»2,...}  to  be  the 


number  of  (a-b)  packets  that  have  arrived  at  origin  a  to  be  transmitted 


th 


over  (a,b)  6  A  to  destination  b  before  the  beginning  of  the  n  slot.  Each 


arrival  process  F  .(•)  :  (a,b)  B  A  is  independent  and  Poisson  with 

V 


rate  f  .  ,  i.e., 


^F*btn+1)  "  rab(n)  "  -  k  >  0,  n  >  1. 


We  define  S  w(n)  n  B  {1,2,...}  to  be  the  number  of  (a  -  b) 

to 


packets  that  have  been  successfully  transmitted  from  origin  a  to 


st 


destination  b  over  (a,b)  before  the  beginning  of  the  (n+1)  slot. 
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Mote  that  a  packet  leaves  the  network  following  its  first  successful  trans¬ 
mission,  so,  f  .  (n)  >  S  .  (n)  for  all  n  >  1,  and  all  (a,b)  8  A. 
an  —  ao  —  . 


A  packet  is  scheduled  for  transmission  in  the  first  slot  following 
its  arrival.  Each  packet  which  is  waiting  for  retransmission  over  link 
(a,b)  is  scheduled  for  the  current  slot  with  probability  0  <  u^  <  1 
independently  of  all  other  packets.  As  a  consequence,  we  admit  the 
possibility  of  more  than  one  packet  being  scheduled  for  transmission  by 
the  same  node  in  the  same  slot.  When  this  happens,  all  packets  involved 
are  treated  as  if  they  have  suffered  collisions  and  the  corresponding 
node  transmits  a  blank  signal  in  that  slot.  Clearly ,  this  system  can  be 
improved  .it  no  cost,  but  we  wish  to  keep  it  this  way  to  facilitate 
the  analysis. 

We  shall  assume  that  acknowledgements  are  available  to  the  trans¬ 
mitting  stations  immediately  following  the  transmissions .  This 
assumption  is  made  for  definiteness  and  can  be  relaxed  without  altering 
our  results. 

We  define  the  scheduling  process  of  packets  over  link  (a,b)  as 
follows: 


Gab(1> 


*ab(1) 


®ab(n+1>  "  Gab(n)  +  Bab(n+1) 


n  6  {1,2,...} 


where  B^tn)  is  the  number  of  packets  scheduled  for  transmission  over 
(a,b)  in  the  n*1*  slot. 


2.3  The  Equilibrium  Hypothesis 

Suppose  the  mean  link  arrival  rates  are  sufficiently  small  so  that. 
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with  high  probability,  packets  are  transmitted  successfully  at  the 
first  transmission  attempt.  Occasionally,  there  will  be  collisions; 
but  if  the  retransmission  probabilities  are  small  enough,  and  if  they 
are  adjusted  with  respect  to  the  prevailing  traffic  density  in  the  network, 
we  expect  resolution  of  contentions  for  channel  vise  and  an  eventual 
return  to  normal  conditions. 

The  above  argument  is  the  basis  of  the  equilibrium  hypothesis, 

which  approximates  the  probability  assignment  for  B^(n)  by 

k  ~ga b 

g  6 

Pte^tn)  -  k}  a  -  (2.1) 

where  the  constant  gaK  is  called  the  mean  packet  scheduling  rate  on  link 
(a,b) . 

The  equilibrium  hypothesis  does  not  hold  for  the  Aloha  model  we  have 
described  in  the  previous  section,  because  if  f^  >  0,  then 

-►  •  asn  +  «  with  probability  1. 

Despite  this  inherent  instability,  the  possibility  of  "stable" 
operation  over  long  periods  of  time  motivates  the  analysis  we  offer 
in  the  next  section. 


2.4  Analysis  of  Slotted  Aloha  Under  the  Equilibrium  Hypothesis 


Let  s^  be  the  probability  that  in  an  arbitrarily  chosen  slot  there 
is  a  successful  packet  transmission  over  (a,b) .  Under  the  equilibrium 
hypothesis,  we  have  that 


or 


*ab 

«ab 


.■*-  TT 

(0,4)*:^ 

•*p{-  £  * 


(2.2) 


(2.3) 


The  summation  in  the  exponent  depends  only  on  node  b  (see  the 

* 

second  remark  following  the  definition  of  C^K  in  Chapter  1).  Therefore, 
we  define 


E. 


(c,d)ec 


*0d 


ab 


and  express  (2.3)  as 


for  all  (a,b)€A  .  (2.4) 


In  equilibrium  we  must  have  s  *  f^  for  all  (a,b)6A.  To  see  why 
this  must  be  so,  suppose  s^h  <  f  for  some  link  (a,b) .  Then,  the 
number  of  packets  waiting  for  retransmission  over  (a,b)  increases 
continually,  thus,  increasing  the  value  of  gah,  which  in  turn  causes 
to  decrease.  As  a  result  equilibrium  cannot  exist. 

On  the  other  hand,  s  >  f  implies  that  the  expected  number  of 

an  ao 

packets  that  are  successfully  transmitted  over  (a,b)  in  a  time  period 
exceeds  the  expected  number  of  new  packets  entering  the  network  for 
transmission  over  (a,b)  in  the  same  time  period.  This  is  clearly 
contrary  co  equilibrium. 

Therefore,  under  equilibrium  conditions  we  must  have 


fab  '  g  -  e 


-*b 


ab 


for  all  (a,b)€A 


(2.5) 


If  we  take  {g^}  as  given,  then  {f^}  is  uniquely  determined  by 


(2.5).  Therefore,  there  are  sets  {f^}  and  {g^}  which  satisfy  (2.5), 
and  the  equilibrium  hypothesis  does  not  lied  lately  lead  to  an  inconsistent 


set  of  equations. 

He  would  like  to  have  necessary  and  sufficient  conditions  on 
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{£  }  which  guarantee  the  existence  of  solutions  to  (2.5) ,  and  which  acre 

ab 

easily  verifiable.  He  do  not  know  of  any  such  necessary  and  sufficient 
conditions  short  of  solving  (2.5). 

However,  a  simple  necessary  condition  on  {f^}  for  the  existence  of 
solutions  to  (2.5)  is  available  through  the  following  observation. 

Suppose  {g^}  is  a  solution  to  (2.5).  Let  us  rewrite  (2.5)  as 
^b“^ab  9ab 

9ab  "  (fab  ®  >  e  ,  for  all  (a,b)eA. 


(2.6) 


The  term  -  g^  in  the  exponent  does  not  explicitly  depend  on  g^, 
thus,  (2.6)  is  an  equation  of  the  form  x  «  aex,  for  which  we  have  three 
cases. 


i)  a  >  e 

ii)  a  ■  e 

iii)  a  <  e 


-1 


-1 


-1 


:  there  is  no  solution. 

:  x  »  1  is  the  only  solution. 

:  there  are  two  distinct  solutions; 

one  greater  than  1,  and  the  other  less  than  1. 


Therefore,  we  conclude  that  (2.6)  has  a  solution  only  if 

Sf’ab  -1 
f»b  0  i* 


or 


"1°  *«b  —  1  “  9ab  +  S> 

Since  g^  >_  f^  for  all  (a,b)€A,  we  have  obtained  the  following 
necessary  condition: 


(2.7) 


(2.8) 


(2.5)  has  a  solution  only  if 

-In  f  .  >1  +  2mJ  V(a,b)8A  s.t.  f  .  >  0  .  (2.9) 

*b  "  (e.diw:^  •*> 
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Notice  that  (2.9)  can  be  satisfied  only  if  f  .  <  e  1  for  all 

an  - 

(a,b)e*. 

In  the  rest  of  this  chapter  we  shall  investigate  a  computational 

algorithm  which  converges  to  a  particular  solution  to  (2.5),  if  one 

exists.  For  this  purpose,  we  find  it  convenient  to  express  our  equations 

-►  -*■ 

in  vector  form.  We  define  g  and  f  to  be  column  L-vectors 
corresponding  to  {g^}  and  {f^} ,  respectively,  and  consider  the 
following  iterative  algorithm. 


g(n+l)  ■  h(g(n) ) 

gd)  -  f 


Z1  8  1  /  2  f  .  •  • 


(2.10) 


where  hab(g(n)>  -  e 


Vn) 


Vn)  * 


E 

(c,d)ec 


*  Wn} 

ab 


A  vector  x  is  called  an  equilibrium  point  of  (2.10)  if  x  ■  h(x). 
The  equilibrium  points  of  (2.10)  correspond  to  t:.e  solutions  of  (2.5). 

We  shall  show  that  (2.10)  converges  to  an  equilibrium  point  if  and 
only  if  (2.5)  has  a  solution. 

Let  x  and  y  be  non-negative  column  L-vectors  such  that  x  ^  y. 

Then,  since  each  row  of  h(<)  is  a  non-decreasinc  function  in  each 
■*  +  +  + 

variable,  h(x)>^  h(y). 

Thus,  in  (2.10),  g(2)  ^  g(l)  -  f  because  g(2)  -  h(f)  and 

.  +  ♦  •*>  > 

g(l)  ■  h(0),  and  f  ^  0.  This  argument  also  shows  that  (2.10)  is 
non-decreasing . 


Let  g  be  any  equilibrium  point  of  (2.10),  if  one  exists. 
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We  have  already  argued  that  g  f.  So,  in  (2.10),  g(l)  £  g.  Now 


4  "4  4-  4  4  4 

suppose,  for  some  n  >_  1,  g(n)  £  g,  then  g(n+l)  =  h(g(n))  h(g)  =  g. 


By  induction,  (2.10)  is  bounded  above. 

Therefore,  (2.10)  converges  (because  it  is  non-decreasing  and 

i 

bounded)  if  and  only  if  (2.10)  has  an  equilibrium  point  (equivalently. 


if  and  only  if  (2.5)  has  a  solution), 


'  4-*  -V 

Suppose  (2.10)  converges  to  g  .  Note  that  g*  is  an  equilibrium  point 


of  (2.10).  The  above  argument  .©hows  that  if  g  is  any  equilibrium  point 

-*•*  / 


of  (2.10),  then  g  >_  g  .  In  pther  words,  if  (2.5)  has  solutions,  there 


exists  a  smallest  equilibrium  point  of  (2.10)  and  (2.10)  converges  to  this 
point. 


-*■* 


In  fact,  if  (2.10)  converges  to  g  ,  then,  starting  from  any  initial 


■4-  4 


point  g(D  in  t^ie  region  0<g  (l)<g* ,  g(n+l)  =  h(g(n))  converges  to  g* .  To 


—  4-4-  -4  4-4" 

see  this,  suppose  0  <  x(l)«  g*  is  arbitrary  and  let  y(l)  ®  0.  Since 


4-  •+  4-  4-*  4-  4-  4-* 

y (2)  »  h(0)  «  f,  y (n)  -*■  g  as  n  +  *  .  But  y(n)  <  x(n)  <  g  for  all  n; 


4-* 


therefore  x(n)  converges  to  g  as  well. 


Suppose  again  that  (2.10)  converges  to  g  .  Let  x(l)  be  such  that 


4*  4-  4-* 

(1  -  e)g  <  x(l)  <  g  where  0  <  e  <  1.  If  e  is  sufficiently  small. 


we  can  approximate  x(n),  defined  recursively  by  x(n)  =  h(x(n-l)) 


(n  >  2) ,  as  follows; 


x(n+l)  ss  g  +  H(x(n)  -  g  ) , 


n  >  1, 


(2.11) 


where  H  is  the  Jacobian  matrix  of  h(«)  evaluated  at  g  . 


By  defining  5(n)  *  g  -  x(n) ,  (2.11)  can  be  expressed  as 


5  (n+1)  ■  H  i  (n) 


n  >  1 


(2.12) 


-i —  .. —  auAttiaB 


4 


Since  5(n)  ■+  $  as  n  -►  »  ,  the  eigenvalues  of  H  must  lie  within 
the  unit  circle  in  the  complex  plane. 

The  elements  of  H  are  as  follows: 


with  this  property,  denoted  by  y  ,  is  called  the  douiinant  eigenvalue 
of  H.  Thus,  (2.5)  has  a  solution  if  and  only  if  y  < 

u  — 

This  argument  is  instructive  and  useful  even  though  y  cannot  be 

u 

computed  easily  from  the  knowledge  of  f  alone.  For  example,  we  can 
-*•*  -*■ 

show  that  g  <  1  is  satisfied  by  making  use  of  some  properties  of  non¬ 


negative  matrices. 

-**  th 

Consider  the  diagonal  matrix  diag(g  )  whose  (a,b)  diagonal  entry 

*  ->.*  -*■* 
is  g  .  >  diag(g  )  is  non-negative  and  H  >  diag(g  ).  The  dominant 
ab  — 

■+* 

eigenvalue  of  diag(g  )  is 


y_  •*  max  {g  ,  :  (a,b)eA}  .  It  follows  that  y_  <  y„ 

G  AD  G  •—  n 

(see,  e.g. ,  [KT]  pages  542-551) . 

Therefore,  (2.10)  converges  to  g  only  if  y_  <  1,  or  in  other 

G  — 

•**  •+. 
words,  if  (2.10)  converges  to  g  ,  then  g  £  1. 
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We  conclude  the  analysis  of  slotted  Aloha  by  obtaining  a  lower  bound 
on  V 

Define  the  matrix  F  as  follows: 


(F) 


ab,cd 


(o.aiec^  ,  fod  >  o 

otherwise 


F  is  obtained  from  G  by  replacing  g^  with  f^  for  each  (a,b)€A. 
Thus,  0  <_ F  ^  H,  and, if  Y_  is  the  dominant  eigenvalue  of  F,  then 

V  l  V 


It  is  clear  that  (2.5)  has  a  solution  only  if  y  <  1. 

C  *“ 

To  obtain  a  lower  bound  on  y  (and  also  on  Y,.)  ,  let  us  define 

F  a 

a^h  and  6^  to  be  the  sum  of  the  elements  in  the  row  and  the  column  of  F, 
both  corresponding  to  link  (a,b) ,  respectively.  That  is. 


-E 

(c,d)ec 


f  . >o 
cd 


£ab  '  f«b  "ab 


where  is  the  number  of  links  in  with  positive  arrival  rates, 

and 


ab 


cd  * 
(c,d)  :  (a,b)6C 


f  .  >  0 
ab 


cd 


f  .  -  0 
ab 
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Let 

Y  *  max{min{a  :  (a,b)€&,  a  >  0}  , 

O  HD  BD 

min  (8^  :  (a,b)eA,  8^  >  0}} 

By  the  general  properties  of  non-negative  matrices  (see,  e.g. , 
page  194  of  [L]),  we  have  that  Yq  5.  YF-  The  computation  of  Y0  i*  easy 
and  (2.5)  has  a  solution  only  if  Yq  1* 

The  results  of  this  section  can  be  summarized  as  follows. 


Proposition  2.1. 

(2.5)  has  a  solution  if  and  only  if  (2.10)  converges.  (2.10) 
converges  only  if 


-1*  fab  i  1  + 


E 

,d)ec 


and  Yp  <,  1. 


V  (a,b)€A,  s.t.  fab>  0, 


2.5  Discussion  of  Results 

As  we  have  remarked  before,  the  Aloha  system  we  have  considered  in 

this  chapter  is  inherently  unstable;  the  main  reason  for  the  analysis 

of  Aloha  was  to  understand  the  conditions  under  which  the  system  could 

be  expected  to  operate  satisfactorily  before  the  throughput  collapsed. 

The  equilibrium  hypothesis  implies  that  stable  operation  is  not 

possible,  even  temporarily,  unless  (2.5)  has  a  solution.  If  (2.5)  has 

a  solution,  stabilisation  of  Aloha  may  be  possible  by  the  introduction 

-► 

of  a  control  mechanism;  in  this  sense  we  can  say  that  f  is  feasible 
under  slotted  Aloha  if  (2.5)  has  a  solution. 


i 


However,  our  results  are  subject  to  the  limitations  of  the 
equilibrium  hypothesis,  which  totally  ignores  the  interdependencies 
among  the  packet  scheduling  rates  of  different  links.  It  would  be 
desirable  to  study  the  f-feasibility  problem  on  a  model  which  at  least 


partially  accounts  for  these  dependencies,  but, mainly  because  of 
tractability  problems,  we  do  not  pursue  this  matter  any  further. 


CHAPTER  III.  TOMA  in  PR *» 

In  this  chapter  we  consider  TDMA  (time-division-multi-access) 
schemes  for  PRNs  After  formulating  TOMA  in  terms  of  the  transedssian 
vectors  of  a  PRN,  we  address  the  feasibility  problems,  which  were 
introduced  in  chapter  1. 

The  main  result  of  this  chapter  is  the  NP-completeness  of  the 
feasibility  problems.  The  reader,  unfamiliar  with  up  concepts,  is 
referred  to  [PS]  or  [GJ]  for  the  essentials  of  algorithmic  complexity. 

3.1  Notation  and  Definitions 

He  change  our  notation  slightly  in  this  and  the  following  chapters 
to  avoid  cumbersome  subscripts. 

As  usual,  G  «  (N,A)  will  represent  a  PRN.  Nodes  will  be  denoted 
by  the  small  case  letters  a,b,c,d  and  e.  L  -  [a|  will  be  the  number 
of  links,  and  the  links  will  be  labelled  by  integers  1  through  L. 

The  letters  i,j,k  and  t  are  reserved  for  denoting  the  links. 

With  the  new  notation,  f^  will  stand  for  the  mean  packet  arrival 
rate  on  link  i.  The  collection  of  all  arrival  rates  will  be  represented 
as  a  set  {f^}  and  alternatively,  as  a  column  L- vector  f ,  the  1th 
component  of  which  equals  f^. 

TOMA  analysis  is  based  on  the  concept  of  collision-free  transmission 
vectors.  To  every  time  instant  we  associate  a  column  L- vector  t  with 
the  following  property:  t^  »  1  if  there  is  a  transmission  over  link  i 
at  that  time  instant;  ■  0  otherwise. 


V 
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A  vector  t  with  the  above  property,  is  called  a  transmission  vector. 

A  collision-free  transmission  vector  t  is  a  transmission  vector 
with  the  following  additional  property:  t^  +  t^  <  1  if  i  6  for 
all  i,j  e  A. 

Since  we  are  interested  only  in  collision-free  transmission  vectors 
in  our  study  of  TDMA,  from  here  on  we  shall  refer  to  collision-free 
transmission  vectors  simply  as  transmission  vectors. 

We  say  that  link  iSA  is  enabled  (or  used)  by  a  transmission 
vector  t  if  »  1.  The  transmission  set  of  a  transmission  vector  £ 
is  defined  to  be  the  set  of  links  that  are  enabled  by  t.  Thus  the 
transmission  set  of  t  is  the  set{i€A:  t^  =  1}  . 

By  convention,  the  all  zeroes  vector  ft  is  a  transmission  vector, 
and  the  null  set  is  a  transmission  set. 

A  maximal  transmission  set  is  a  transmission  set  which  is  not 
properly  contained  in  another  transmission  set.  A  maximal  trans¬ 
mission  vector  is  a  transmission  vector  whose  transmission  set  is 
maximal. 

The  total  number  of  all  transmission  vectors,  denoted  by  K,  is 
obviously  bounded  by  2L.  We  let  t^,^,...,^  be  an  ordering  of  all 
transmission  vectors. 

The  transmission  matrix  T  is  defined  to  be  an  LxK  matrix  whose 
1th  column  is  for  i  -  1,...,K. 

In  the  L-dimensional  read  space  RL,  the  transmission  vectors 
can  be  thought  of  as  points.  The  convex  hull  of  all  those  points, 
denoted  by  CH(T),  is  a  closed  bounded  poly tope  which  lies  in  the  non¬ 
negative  orthant  of  RL. 


3.2  TDMA  in  PBHs 

M  are  now  in  a  position  to  introduce  the  TDMA  scheme.  For 
simplicity  of  exposition,  in  this  chapter  we  consider  variable  slot 
length  TOMA  schemes.  In  the  variable  slot  length  case,  each  user  has 

a  sequence  of  positive  numbers.  The  sequences  are  identical.  The 

* 

terms  of  the  sequence  represent  the  duration  of  slots.  All  users  are 
synchronised  so  that  the  slots  as  perceived  by  different  users  always 
start  and  end  simultaneously. 

In  TOMA,  we  associate  a  transmission  vector  with  every  slot.  The 
sequence  of  transmission  vectors  associated  with  the  slots  is  assumed 
to  be  the  common  knowledge  of  all  users.  In  our  study,  we  do  not 
distinguish  the  order  in  which  the  transmission  vectors  are  used. 

If  in  some  period  of  time  the  same  transmission  vector  is  used  more 
than  once,  we  caui  combine  all  the  slots  in  which  it  is  used  into  a 
single  larger  slot  whose  duration  is  the  sun  of  the  durations  of 
the  individual  slots. 

Since  we  are  studying  the  feasibility  problems  in  a  static  setting 

we  restrict  our  attention  to  periodic  sequences  of  transmission  vectors 

Me  let  x^  denote  the  duration  of  time  the  transmission  vector  is 

used  in  one  period.  The  column  K-vector  x  will  denote  the  vector 
til 

whose  i  row  is  x^.  We  put  a  normalization  constraint  on  x  by 
demanding  that  t  •  x  <  1. 

After  this  introduction,  we  give  a  formal  definition  of  TDMA. 
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Def inltion .  Let  G  ■  (H,A)  be  a  PRN  and  let  T  be  the  transmission 

matrix  of  G.  A  TDMA  scheme  (for  G)  is  a  system 
<T ,x> ,  where  x  is  a  non-negative  column  K-vector 
such  that  t  •  x  <  1. 

In  a  PKN,  all  links  have  the  same  capacity.  We  normalise  link 
rates  with  respect  to  this  capacity  and  require  that  0  <_  ?  <_  Ir . 

Definition.  Let  G  -  (N,A)  be  a  PKN  and  let  1  be  a  link  rate  vector 
for  G.  1  is  called  feasible  under  TDMA  if  there 

^  —V 

exists  a  TDMA  scheme  <T,x>  such  that  Tx  >  f , 

We  will  now  state  a  series  of  straightforward  facts  about  TDMA 
schemes.  Pacts  3.1  -  3.5  are  stated  in  reference  to  a  fixed  but 
arbitrary  PRN  G-(N,A). 

FACT  3.1. 

?  is  feasible  under  TDMA  iff  u  <  1  where  u  is  defined  as  follows: 
LPl: 

u  -  min  t  •  x 
s.t.  TX  >. 

FACT  3.3 » 

If  is  feasible  under  TDMA,  then  there  exists  a  TDMA  scheme 
<T,x>  such  that  Tx  >  ?  and 

|{i  »  xL  >  0>j  <  L. 
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Proof.  If  f  is  feasible  under  TDMA,  then  there  exists  a  basic 

feasible  solution  which  optimises  LPl  (the  linear  program 
defined  in  the  statement  of  Fact  3.1). 

FACT  3.3. 

Let  T  be  the  set  of  all  transmission  sets  of  G  =  rN,A) .  The 
•  system  (A,x)  is  an  independence  system. 

Proof .  A  system  (E,S)  is  an  independence  system  if  S  is  a  set  of  sub¬ 
sets  of  E  with  the  following  property:  E^es  and  implies 

E2es. 

Since  a  subset  of  a  transmission  set  is  a  transmission  set,  (A,T) 
is  an  independence  system. 

FACT  3.4. 

Let  u  be  as  given  in  Fact  3.1,  let  v  be  as  follows: 

LP2: 

v  ■  min  1  •  x 
s.t.  Tx  »  f 
x>0 

Then,  u  ■  v. 

Proof.  Fact  3.4  follows  from  the  independence  system  property  of  (A,T) . 
FACT  3.5. 

1  is  feasible  under  TDMA  iff  f  6  CK(T). 
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Proof.  £  6  CH(T)  iff  f  can.  be  expressed  as  a  convex  combination  of  the 
extreme  points  of  CH(T),  which  are,  by  definition,  the  transmission 
vectors  of  G. 

FACT  3.6. 

-> 

If  f  is  feasible  under  TCMA,  then  there  exists  a  TOMA  scheme 
-*■ 

■  <T,x>  such  that 

(i)  T*>f, 

(ii)  j  e{i  :  >  0}  t^  is  maximal, 

(iii)  |{i  :  xL  >  0} J  <  L. 

FACT  3.7. 

Given  a  polynomial  p(  •),  we  can  always  find  a  PRN  with  the 
property  that  K  >  p(L) . 

Proof.  An  example  suffices.  Consider  the  PRN  in  Figure  3.1  with 

(1  +  ~)  nodes  and  L  links  (I»  is  even) .  There  is  a  trans¬ 
mission  set  with  L/4  elements.  Therefore,  the  total  number 
of  transmission  vectors  is  at  least 


A  Chain  PBN 
Figure  3.1 


4 
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As  a  result,  linear  programming  as  applied  to  LPl,  or  LP 2  does 
not  guarantee  an  efficient  solution  to  the  f-feasibility  problem. 

Can  we  reduce  the  size  of  LPl  (or  LP2)  by  discarding  non-maximal 
transmission  vectors,  so  that  linear  programming  can  be  used  to  solve 
the  f-feasibility  problem  efficiently?  The  following  fact  suggests 
that  we  cannot . 

FACT  3.8. 

Given  a  polynomial  p(*),  we  can  always  iind  a  PRN  with  the  property 
that  X"*  >  p  (L) ,  where  x'  is  the  number  of  all  maximal  transmission 
vectors. 

Proof.  Consider  a  chain  PRN  G  ,  such  as  the  one  in  Figure  3.1  with 

—  m 

3m+2  nodes  and  L  ■  6m+2  links  for  some  integer  m  >_  1.  By 

induction  we  will  show  that  G  has  at  least  2*3m  maximal  trans- 

m 

mission  vectors. 

G  .  is  constructed  by  adjoining  three  nodes  to  G  ,  as  in  Figure  3.2- 
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Let  S  be  a  maximal  transmission  set  of  G  .  There  axe  four  cases: 

SI 

1)  (3m+l,  3m+2)  e  S 

In  this  case  SU{(3m+4,  3m+5)}  ,  S U (  ( 3m+4 ,  3m+3)  }  , 

and  SU((3m+5,  3m+4) }  are  maximal  transmission  sets  of  G  ... 

StX 

2)  (3m+2,  3m+l)  6  S 

In  this  case  SU{(3m+3,  3m+4)}  ,  S(j{(3m+4,  3m+5) }  ,  and 

sU{(3m+5,  3m+4) }  are  maximal  transmission  sets  of  G  . 

m+i 

3)  (3m,  3m+l)  e  S 

In  this  case  SU{(3m+3,  3m+2),  (3m+4,  3m+5)},  SU(  (3m+3,3m+4) }, 

and  S  U  {(3m+5,  3m+4) }  are  maximal  transmission  sets  of  G 

—  m+1 

4)  (3m+l,  3m)  e  S 

In  this  case  St>{(3m+2,  3m+3) ,  (3m+5,  3m+4)}  ,  SU{(3m+3,  3m+4)} 

SU{(3m+4,  3m+3)}  ,  and  SU{(3m+4,  3m+5) }  are  maximal 

transmission  sets  of  G  ... 

m+1 

Since  every  maximal  transmission  set  of  G  must  fall  into  one  and  only 

81 

one  of  the  four  cases  listed  above,  we  have  shown  that  G ■  .  has  at  least 

three  times  as  many  maximal  transmission  sets  as  G  . 

m 

For  m  ■  1,  by  enumeration,  we  see  that  there  are  at  least  2*31  maximal 

transmission  sets.  Therefore,  we  conclude  that  G  has  at  least  2*3m 

m 

maximal  transmission  sets. 

In  order  to  discuss  the  complexity  of  the  ^-feasibility  problem 
under  TDMA,  we  need  to  define  a  size  for  the  problem.  This  forces  us  to 
make  the  following  assusqption.  The  link  arrival  rates  are  rational 
numbers  such  that  there  exists  a  fixed  (but  arbitrary)  positive  integer  P 
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with  the  following  property:  P*f^  is  integer  for  all  i€A,  It  is 
important  to  note  that  P  stays  fixed  over  all  instances  of  the  problem, 
but  otherwise,  P  is  arbitrary. 

Definition. 

FF  (TDMA  -  f  -  Feasibility  Problem) 

Instance:  A  PRN  G  =  (N,A),  a  link  arrival  rate  vector  f. 

♦ 

Question:  Is  f  feasible  under  TDMA? 

We  define  the  size  of  an  instance  of  FF  to  be  L. 

We  now  state  the  major  result  of  this  section  as  a  theorem,  the 
proof  of  which  is  deferred  to  the  next  chapter. 

Theorem  3.1. 

FF  is  NP-complete . 


The  r-Feasibility  Problem 


Lei  r  be  an  o-d  arrival  rate  vector  of  a  PRN,  G  *  (N,A).  We  assume 
that  r  is  rational?  that  is,  there  exists  a  fixed,  but  arbitrary, 
integer  P  such  that  P  r  is  integer. 


Definition. 

RF  (r-Feasibility  Problem) 

Instance:  A  PRN  G  *  (N , A) ,  an  i-d  arrival  rate  vector  r. 

Question:  Is  r  feasible,  that  *«,  is  it  true  that  r  e  C(G)? 

The  above  problem  is  the  same  as  the  problem  defined  ii:  section  3 
of  Chapter  1  except  for  the  assumption  that  r  is  a  rational  vector. 


We  shall  discuss  the  significance  of  this  assumption  in  the  next  chapter. 
We  close  this  section  by  stating  a  theorem  which  will  be  proved  in 
Chapter  IV. 

Theorem  3.2. 

RF  is  NP-complete. 

3.4  Discussion  of  Results 

The  fact  that  such  fundamental  problems  as  FF  and  RF  are  NP-complete 
is  discouraging,  because  the  lack  of  adaptability  of  the  network  protocols 
to  unexpected  changes  in  network  topology  or  in  data  traffic  makes  proper 
operation  of  the  PRN  uncertain. 


CHAPTER  IV.  NP  COMPLETENESS  PROOFS  AND  AN  ALTERNATIVE  FORMULATION  OF  TOMA 

In  this  chapter  we  shall  prove  the  theorems  we  stated  in  Chapter  III, 
and  give  a  formulation  of  TDMA  for  the  case  where  all  slots  sure  the  same 
in  length. 

The  transformation  algorithms  we  use  in  this  chapter  and  the  next 
require  constructions  of  certain  sets  from  other  set3.  The  set  notation 
that  we  use  in  these  algorithms  is  not  a  standard  one,  so,  we  should 
explain  the  notation  first. 

Let  V  be  a  set,  say  V  «  {x,y7z},  then  V1*  denotes  the  set 
{x",y",  z'}.  In  general,  we  write  V'  =  {a "j  a  e  v}  .  Likewise,  V" 
is  the  set  {a"*  :  a  6  V}  . 

Often,  we  use  indices  instead  of  primes.  If,  for  example, 

X  =  {a,  b},  then  Y1  =>  {a^  b  }  ,  Y2  “  {a2,  b2>  ,  and  so  on.  In  general, 

we  write  Y.  =  {a. :  a  6  Y}  . 

11  J 

This  notation  applies  to  sets  of  ordered  pairs  as  well.  For 

\ 

example,  let  E  =  {(a,b),  (c,d),  (a,d)},  then  tl^e  set  :  (x,y)  6  E} 

written  by  its  elements  is  {(a^  b^ ,  > .  Similarly, 

{<xi»y^)  *  <x,y)  6  E}  ■  {(ai,b^),  (cif  d^) ,  (a^,  d^)}  . 

Finally,  let  us  note  the  following  type  of  equivalences: 

{a,  a"  :  a  8  V>  ■  VUV'  , 

{ (a,b') ,  (a',b)  ,  (a,  b)  6  E>  - 

{ {a,b")  :  (a,  b)  e  e}  U  {(a',b)  :  (a,b)  6  E>  . 
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In  this  chapter,  H  *»  (V,E)  denotes  an  undirected,  finite,  and 
simple  graph. 

4.1  NP-Comp leteness  of  FF 

The  NP-completeness  proof  for  FF  will  be  given  in  two  steps.  First, 
we  shall  consider  a  related  problem  and  prove  that  it  is  NP-complete . 

Then,  it  will  be  shown  that  this  problem  is  polynomially  transformable 
to  FF. 


4.1.1  NP-Completeness  of  the  Maximum  Transmission  Set  Problem 
Definition.  MTS  (The  Maximum  Transmission  Set  Problem) 

Instance:  PRN  G  =  (N,A)  set  B  C  A,  positive  integer  k  <  |  b}. 
Question:  Does  there  exist  a  transmission  set  S  of  G  such 
that  |  S  O  B  |  «  k? 

We  will  prove  the  NP-completeness  of  MTS  by  polynomially  trans¬ 
forming  the  following  NP-complete  problem  to  MTS. 

Definition.  MC  (The  Maximum  Clique  Problem) 

Instance:  Graph  H  =*  (V,E),  positive  integer  k  £  |v|  . 

Question:  Does  G  contain  a  clique  of  size  k  ? 

(A  clique  of  H  ■  (V,E)  is  a  set  Q  £V  such  that  (a,b)  6  E  for  every 
pair  of  distinct  nodes  a,  b  e  Q) . 

Fact  4.1.  MC  is  NP-complete. 

Proof.  See,  e.g.,  page  360  of  IPS) . 


MC  is  polynomial ly  transformable  to  MTS  by  using  the  following 


algorithm. 


Algorithm  Al. 

Input:  Graph  H  «  (V,E)  ' 

Output:  PEN  G  -  (N,A) ,  set  A  where 
N  -  V'kJV" 

A  -  { (&' ,b") ,  aev,  bev,  (a,b)  0  E} 

B  -  {  (a' ,  *")  :  aev} 

We  illustrate  the  transformation  by  an  example. 


Example.  Let  H  »  (V,E)  be  as  in  Figure  4.1. a. 

G  -  (N,A) ,  the  output  of  Al,  is  shown  in  Figure  4.1.b 

l 

B-{(a%a~),  <b',b~),  (c'.c')}  . 

G  -  (N,A) 


H  -  (V,E) 


Figure  4,1 
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Fact  4.2.  Let  H  ■  ,{V,E)  and  G  =  (N,A)  be  the  input  and  the  output  of 

Al,  respectively.  H  has  a  clique  of  size  k  iff  there  exists  a 
transmission  set  S  of  G  such  that  |sH  b|  ■  k. 

Proof.  (i)  Suppose  Q  is  a  clique  of  H  with  |q|  »  k. 

Then,  S  **  {(a  ',a  ")  :  a6Q)  is  a  transmission  set  of  G  and  JsH  B  |  «  k. 

(ii)  Suppose  S  is  a  transmission  set  of  G  with  the  property 

that  |  S  O  b|  -  k.  Then,  Q  -  {a  e  V  s  (a',a~)  e  (SO  B)>  is  a  clique  of 
H  and  |q|  *  k. 


Fact  4.3.  MTS  is  UP  -  complete. 

Proof.  It  is  clear  that  the  set  S  is  a  concise ' certificate  that  can  be 
checked  in  polynomial  time  for  validity  for  a  YES  instance  of  MTS. 

By  Facts  4.1  -  4.2,  we  know  that  MTS  is  as  hard  as  any  problem  in  NP. 
Therefore,  MTS  is  NP-complete . 


4.1.2  MTS  is  Polynomially  Transformable  to  FF. 
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{(a^) 


a.e  N. ,  b.e  N. 
i  i  x  1 


ai  *  V 


If 


i  *  j 


A  .  «  {(a.  ,b  )  »  a.e  N.,  b.e  N.  ,  a  ■  b  or  (a,b)6  a}  i  -  l,...,m 
13  13  1  3  3  j  -  l,...,m 

The  best  way  to  understand  the  structure  of  Gn  is  to  consider  an 
example. 

2 

Example.  Let  G  be  as  in  figure  4.2a.  G  is  the  graph  in  Figure  4.2.b. 
Verify  the  following  on  Figure  4.2. 


1)  (N^,  A^)  and  (N2 'a22^  are  complete  directed  graphs. 

2)  (b,a)  e  A  does  not  conflict  with  (e,d)  e  A  in  G;  (b1#a^)  €  A^ 

2 

does  not  conflict  with  (c^  ,d2>  in  G  .  More  generally,  if 
(x,y)  e  A  does  not  conflict  with  (ufv)  6  A,  then  (x^.y^  6  A^ 
does  not  conflict  with  <U2»V2^  6  A22J  is  4130  true  that 
(x2,y2)  does  not  conflict  with  (u^v^. 

With  the  above  observations,  we  can  prove  the  following  general 
result  about  the  powers  of  G. 

Fact  4.4.  Let  G  **  (N,A)  and  B  C  A  be  given,  let  G*  be  the  k**  power  of 
k  * 

G,  and  let  B  -  (J  B  where 

i-1  1 

B±  -  i (ai,bi)  i  (a,b)  e  B> 

Then,  the  following  statements  are  equivalent. 


(i)  There  exists  a  transmission  set  S  of  G  such  that 
|SD  b|  -  k. 
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(A  link  without  an  arrow  stands  £or  two  oppositely 
links  in  this  figure.) 


(b)  G2  -  (N2,  A2) 


Figure  4.2 
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k  k 

(ii)  There  exists  a  transmission  set  S  of  G  such  that 

|skn  3k|  -  k. 

Proof .  Suppose  (i)  is  true .  Let 

S  OB  ■  Let  ^ij  ^  tke  copy  of  ^i  *n 

A^#  that  is,  if  -  (a,b),  then  4.^..  »  {a^,b^)  for  j  ■  l,...,k. 

)c  1  Jc 

Then,  S  *»  { £  :  i  -  l,...,k}  is  a  transmission  set  of  G  , 

)t 

because,  for  any  two  distinct  links  i  e  S  ,  if  conflicts 

with  £^j»  then  ^  fi  S  conflicts  with  e  Si  thus,  S  cannot  be  a 
transmission  set  contrary  to  the  hypothesis. 

Sk  contains  one  link  In  (i  «  l,...,k);  hence.  |skn  Bkj  *  |sk|  «  k 

Suppose  (ii)  is  true.  Then  we  note  that  Sk  contains  one  link  in 
B^  (  i  ■  l,...,k),  because  any  two  distinct  links  in  Bi  (i  *  1,.... ,k) 
conflict  with  each  other.  In  fact,  S  cannot  contain  more  than  one 
link  in  A^  (i  -  l,...,k).  Thus, 

S  «  {(a,b)  :  (a^,b^)  6  S  for  some  i  e  {l,...,k}}  has  exactly 

k  elements.  Moreover,  SC  B  and  S  is  a  transmission  set  of  G. 

This  completes  the  proof. 

MTS  is  polyncmially  transformable  to  FF  by  using  the  following 
algorithm. 

Algorithm  A2. 

*«*  ^ 

(The  graph  G  -  (N,A)  which  appears  at  the  output  of  this  algorithm 
has  a  special  structure  which  can  be  visualized  easily  by  first  skimming 
over  the  following  definitions,  skipping  the  part  about  f  ,  studying 
the  example,  and  then  carefully  reading  the  algorithm  again.) 
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zij  "  ^(Pi'aj)»  *•  Pi  6  aj  6  Nj>  5  *  1  i 

i  ■  1| « •  •  |k  s  3  m  1/  •  •  *  •  /k 


— 

f 

xy 


6 

(  |b  |  -  1)  6 

0 


(x,y)  e  B* 
(x,y)  6  D 
otherwise 


where 


5 


1 

1  +  k  (|Bj  -  1) 


B 


u 


Bi  "  ^(ax'bi)  :  ai  6  Ni'  bi  6  Ni'  (a'b)  6  B*  1  “  1'***'k 

D  *  :  i  «  1, . . .  ,k) 

Example.  Consider  an  instance  of  MTS  <G,B,k>  where  G  *  (N,A)  is  a 

PRN,  B  G  A  is  a  subset  of  links,  and  k  <  |b[  is  a  positive  integer. 
Suppose  |b(  ■ 4  and  k  »  3. 

With  <G,B,k>  as  the  input,  the  output  of  A2  is  illustrated  in 
Figure  4.3. 

In  Figure  4. 3. a,  we  see  the  general  structure  of  G  -  (N,A) .  The 

%  3 

part  of  G  that  is  anclosad  in  tha  da  shad  ractangla  is  the  saa»  as  G  * 

v 

is  shown  in  Figure  4.3.bi  we  note  that  (P,Q,y  )  is  a  complete 

v 

bipartite  directed  graph,  that  is,  <x,y)  6  Y  iff  either  x  6P  and  y  eQ , 


or  x  6  Q  and  y  €P .  in  other  words,  each  node  in  P  is  connected  to  each 
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node  in  Q,  and  vice  versa.  Thus,  a  transmission  set  S  of  G  can 

]c 

have  at  most  one  link  in  common  with  Y  . 

The  cross-hatched  area  in  Figure  4. 3.c  represents  the  links  in 
Z^(if*5);  we  note  that  {{p^},  Z„)  is  also  a  complete  bipartite 

directed  graph.  Hence,  a  transmission  by  node  interferes  with  a 
simultaneous  transmission  over  any  link  in  A^  (j^i).  So,  if  a 

—  **  <v  ** 

transmission  set  s  of  G  contains  a  link  (p^ ,x)  6  A,  then  S  cannot 
contain  any  link  in  A^  (j?*i). 

Finally,  in  Figure  4.3.d,  we  see  the  links  in  D.  Bach  of  the 
links  in  D  is  assigned  a  mean  link  arrival  rate  of  ( [ B )  -  1)  6  *  0.3 
units.  (Here  we  do  not  use  packets/slot  as  the  unit  of  arrival  rates, 
because  slot,  as  a  time  unit,  does  not  have  a  meaning  in  a  variable 
slot  length  TDMA  scheme.  Instead,  we  interpret  the  arrival  rates  as 
service  requirements;  thus,  a  mean  link  arrival  rate  of,  say,  0.1  units 
implies  that  the  link  in  question  must  be  enabled  0.1  fraction  of  time 
to  satisfy  its  service  requirement;  sometimes,  we  express  the  same 
fact  by  saying  that  the  link  has  a  service  requirement  of  0.1  units.) 

The  following  result  is  based  on  the  properties  of  transmission 
sets  of  G  which  have  been  pointed  out  in  the  above  example. 

Fact  4.5. 

Let  the  instance  <G,B,k  >  of  MTS  be  the  input  of  A2,  and  let 
<G,f>  be  the  corresponding  output.  Then,  the  following  statements  are 
equivalent. 

(i)  There  exists  *.  transmission  set  S  of  G  such  that  |s  H b|  *  k. 

(ii)  {f  t  (x,y)  €  A)  is  feasible  under  TDMA. 
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Proof.  There  are  k  |b  J  links  in  each  with  5  units  of  service 

requirement,  i.e.,  each  link  in  B  must  be  enabled  6  fraction  of  time 
on  the  average;  otherwise,  packets  start  piling  up.  There  are  k  links 
in  D  each  with  a  service  requirement  of  (  (b  |  -  1) 6  units.  So,  the 
sum  of  the  service  requirements  on  all  links  of  G  is  k5  +  2k  ([b|  -1)6 
units  (or  k5  +  2(1  -  6) units). 

Now,  suppose  (i)  is  true.  Then,  there  exists  a  transmission  set 
S  of  G  such  that  |sH  b  |  =  k.  s  can  be  used  6  fraction  of  time  to 

jc 

satisfy  the  service  requirements  of  k  links  in  B  .  So,  if,  in  the 
remaining  (1  -  5)  fraction  of  time,  we  can  complete  the  service  of 
k ( j  Em  -  1)  links  in  B  each  with  a  requirement  of  6  units,  and  k  links 
in  D  each  with  a  requirement  of  (|b|  -  1)6  units,  then  statement  (ii) 
will  be  true.  We  shall  show  that  this  is  indeed  the  case. 

We  know  that  S  has  one  link  in  each  (i=l,...,k).  So,  let 
us  write 

B^  -  S  = ’•  j  =  l,...,k-l}  i  ■  l,...,k.  The  set 

S, .  ■'  {(p. ,q.).  A..}  is  a  transmission  set  of  G  for  each  A.,  e  (B.  -  S) , 

i}  31  ji  i 

i  *  1, . . . ,k.  ■ 

It  is  easy  to  see  that  the  TDMA  scheme  which  uses  S  and  S^ 

(i  *  l,...,k  1  j  »  l,...,k-l)  each  for  6  fraction  of  time  satisfies 
the  service  requirements  on  all  links  of  G  just  on  time. 

Suppose  (ii)  is  true.  We  note  that  exactly  k(|B|  -  1)6  (or  (1-6)) 
fraction  of  time  must  be  spent  in  satisfying  the  service  requirements 
of  the  links  in  D.  In  the  remaining  6  fraction  of  time,  the  TDMA  scheme 
under  which  {f  }  is  feasible  must  satisfy  a  total  of  at  least  k6  units 
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o£  service  requirements  on  the  average.  (This  last  sentence  follows  by 
noting  that  the  TDMA  scheme  cannot  use  transmission  sets  of  size  larger 
them  2  in  the  (1-6)  fraction  of  time  that  it  spends  while  serving  the 

jc 

links  in  D) .  Therefore,  the  average  number  of  links  in  B  that  must  be 
enabled  by  the  transmission  sets  used  in  the  remaining  6  fraction  of 
time  is  at  least  k.  This  implies  that  there  exists  a  transmission  set 
S  of  G  such  that  js  OB'  |  =  k.  But,  the  existence  of  S  implies  that 
there  exists  a  transmission  set  S  of  G  such  that  |s  D  b|  =  k.  This 
completes  the  proof. 

To  prove  that  FF  is  NP-complete,  we  need  to  show  that  every  YES 

instance  of  FF  has  a  concise  certificate  which  can  be  checked  in  time 

bounded  polynomially  in  the  size  of  the  instance. 

Let  1  be  feasible  under  TDMA.  There  exists  linearly  independent 

transmission  vectors  t  ,t  , ...,t  such  that,  for  some  x  >  (J,  Bx  =  r 

X  2  L  — 

th  "*■ 

and  1  •  x  1  where  B  is  the  LxL  matrix  whose  i  column  is  t^  (i**l, . . .  ,L) . 

We  propose  B  and  x  as  a  concise  certificate  for  FF.  First,  let  us 
see  if  x  can  be  encoded  concisely. 

-►  -i 

We  have  x  =  B  f. 


Let  (B  )^  be  the  i-j^1  element  of  B  \ 


that  is. 


) 


ij 


adi(Bu> 

det  B 


Since  adj(B^)  is  the  determinant  of  an  (L-l)  x(L-l)  0-1  matrix, 
|adj(B^)  |  '<  (L-2)  !  Likewise,  |det  b|  <  (L-l)!.  Therefore,  each  of 
|det(B) |  and  |adj(B^) |  can  be  encoded  in  fewer  than  L  bits.  Since  f 
is  assumed  to  be  rational,  we  conclude  that  x  can  be  encoded  in  space 


1 


bounded  polynomially  in  L. 

The  columns  of  B  can  be  checked  in  polynomial  time  for  being 
transmission  vectors.  The  operation  Bx  can  also  be  carried  out  in 
polynomial  time.  Therefore,  if  f  is  feasible  under  TDMA,  then  there 
exists  a  concise  certificate.  This  proves  that  FF  belongs  to  NP. 


Notice  that  we  have  not  used  the  assumption  that  f  is  rational  in 
proving  the  fact  that  MTS  is  polynomially  transformable  to  FF.  This 
assumption  was  only  necessary  to  show  that  FF  belongs  to  NP. 


4.2.  NP  -  Completeness  of  RF 


We  define  the  total  o-d  rate  associated  with  node  x  of  PEN 
G  *  (N,A)  as  follows: 


X  XV 

y€N 


Rx  can  be  interpreted  as  the  fraction  of  time  node  x  is  receiving 

or  transmitting  packets  which  either  originate  from  x  or  leave  the 

network  at  x.  If  r  is  feasible,  then  R  <  1  for  all  x€N. 

x  — 

The  following  algorithm  polynomially  transforms  FF  to  RF. 


Algorithm  A3. 


Input:  PRN  G  -  (N,A) ,  set  {f  :  (a,b)  6  A} 

Output:  PRN  G  -  (N,A) ,  set  {r  :  xSN,  yeN> 

xy 


where  N  {a, a", a"  :  a  e  N} 


A  -  A  \j  Ax  U  A2 

Ax-  { (a", a"')  :  a  6  N} 
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A2“{(afa'^>,  (a", a),  (a",aj,'»  (a,a') ,  {&",&')  i  a  e  N> 


rab  jab 


r  >  «  1  -  T  (f.+f.  ) 

a  a  ab  ba 

~  D 

r  =0 

xy 


(a,b)  e  A 

(a',a")  e  Ax 
(x,y)  4  A  or  (x,y)  e  A2 


Note  that  r--.-  =  l-  R'.  Without  loss  of  generality,  we  shall  assume 

S  3  cl 

that  0<  R  <  1  for  all  x  6  N.  When  {f  .  :  (a,b)  e  a}  is  feasible  under 
—  x  —  (  ab 

TOMA,  this  assumption  will  always  be  true. 

We  illustrate  the  algorithm  by  an  example. 


Example.  Let  G  =  (N,A)  be  the  PRN  in  Figure  4. 4. a. 
**■’  «*•  ^ 

G  «  (N,A)  is  shown  in  Figure  4.4.b. 


(a)  G  -  (N,A) 


rab 

m 

0.2 

r* 

rba 

rbc 

m 

0.1 

r* 

rcb 

*+ 

r 

B 

0.5 

r  -  -- 

ca 

a  a 

.* 

r 

m 

0.4 

r  -  -- 

b  b 

c  c 

All  other  o-d  rates  are  zero 


Figure  4.4 


A 
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Fact  4.6.  Let  _.a  PRN  G  =  (N,A) ,  and  the  set  of  mean  link  arrival  rates 

•a 

{f  .  s  (a,b)  6  A}  be  the  input  of  A3,  and  let  G  =  (N,A)  and 
ab 

{r  :  x  8  N,  y  S  N}  be  the  corresponding  output.  The  following  statements 

xy 

are  equivalent. 

(i)  {f  .  :  (a,b)  e  A}  is  feasible  under  TDMA 

ab 

(ii)  .  {r  :  x  e  N,  y  e  N}  is  feasible, 

xy 

Proof.  Suppose  there  exists  a  TDMA  scheme  for  G  under  which 

{f^k  s  (a,b)  e  A}*  is  feasible.  We  shall  show  that  {r^  s  x  &  N,  y  e  N} 

is  feasible  in  G.  Since  r  =  0  when  (x,y)  0  A,  each  packet  in  the 

xy 

network  cam  be  transmitted  directly  from  its  origin  to  its  destination. 

Let  this  be  our  rule  for  the  assignment  of  paths  to  packets;  that  is, 

every  (x  -  y)  packet  is  transmitted  directly  over  link  (x,y)  for 
—  — 
all  (x,y)  e  A.  With  this  path  assignment,  we  have  f^  ■  r^  for  all 
«• 

(x,y)  e  A,  where  f  is  the  mean  arrival  rare  of  link  (x,y) . 
xy 

If  S  is  a  transmission  set  of  G  such  that 

SO  { (x,y)  e  A  :  x  =»  a  or  y  =  a}  -  0  , 

then  S  u{(a',a"')}  is  a  transmission  set  of  G.  Thus,  given  a  TDMA 
scheme  for  G  under  which  {f  .  :  (a,b)  e  A}  is  feasible,  we  can 
augment  the  transmission  sets  used  by  this  TDMA  scheme  to  obtain  a 
TDMA  scheme  for  G  under  which  {f  s  (x,y)  e  A}  is  feasible. 

xy 

Suppose  (ii)  is  true.  Consider  nodes  a,  a',  a"  6  N  .  Node  a " 

must  be  transmitting  at  least  x  ,  ,,  fraction  of  time. and  node  a  cannot 

a  a 

transmit  or  receive  from  nodes  other  than  a "  during  this  time.  Node  a 

must  spend  at  least  1-r  ..  ,,  m  R  fraction  of  time  receiving  from  and 

a  a  a 
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transmitting  to  nodes  other  than  a'  and  &"  .  Thus,  all  (a"  -  a")  traffic 

must  be  transmitted  directly  over  (a",  a") .  But,  if  all  (a"  -  a") 

traffic  is  sent  directly  over  {.a',  a"),  node  a  cannot  handle  any 

packets  other  than  those  which  either  originate  at  node  a  or  leave  the 

network  at  node  a.  Therefore,  one  must  be  able  to  send  each  packet 

^  -  ** 

directly  from  its  origin  to  its  destination  if  {r  :  x  e  N,  y  G  N } 

xy 

is  feasible.  This  completes  the  proof. 

A  concise  certificate  for  a  YES  instance  of  RF  is  a  set  of  routing 
variables  and  a  TDMA  scheme  under  which  the  induced  set  of  mean  link 
arrival  rates  are  feasible. 

Thus,  we  have  proved  that  RF  is  NP-complete. 

4.3  Fixed  Slot  Length  TDMA  Schemes 

Until  now,  we  have  considered  variable  slot  length  TDMA  schemes. 
Sometimes,  it  may  be  desirable  to  have  a  scheme  with  all  slots  fixed 
in  length.  In  this  section,  we  shall  show  that  this  can  be  done  without 
loss  of  capacity,  in  the  sense  that,  if  a  set  of  mean  link  arrival  rates 
are  feasible  under  TDMA,  then  there  exists  a  fixed  slot  length  (FSL; 

TDMA  scheme  under  which  the  same  set  of  mean  link  arrival  rates  are 
feasible.  The  converse  of  this  statement  is  also  true. 

«► 

Definition.  Let  G  ■  (N,A)  be  a  PRN,  and  let  f  be  a  mean  link 
arrival  rate  vector .  A  FSL  -  TDMA  scheme  with  period  J  is  a  system 
<T,  u,J>  where  T  is  a  transmission  matrix  of  G,  u  is  a  column 
k  -  vector  with  non-negative  integer  elements  and  J  is  integer  such 
that  J  >  1  •  u.  We  say  f  is  feasible  under  FSL-TDMA  if  there  exists 
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a  FSL-TDMA  scheme  <T,u, J>  such  that  Tu  >  J  f 


Fact  4.7.  If  for  some  P,  P  •  f  has  integer  elements,  and  if  f 

is  feasible  under  TDMA,  then  there  exists  a  FSL-TDMA  scheme  <T,u,J> 

+  • 
such  that  T  u  >  J  f  and  3  <  P  •  (L  -  1)  .  , 


-T* 

Proof.  If  f  is  feasible  under  TDMA.,  then,  as  we  have  shown  in  the 

t 

NP-completeness  proof  of  FF,  there  exists  a  TDMA  scheme<T,x>  such  that 

-v 

Tx  **  f  and  P*  Mx  has  integer  elements,  where  M  is  the  determinant 
of  an  LxL  0-1  matrix.  Since  the  determinant  of  an  LxL  0-1  matrix 

must  be  less  than  (L-l) !  ,  we  have  that  M  <  (L-l) !  .  Let  u  •>  P*  M  x 

^  ^  ^ 
and  J  =  P  •  M.  Then,  Tu  =  P  •  M  »Tx*»P*Mf**J  •  f  and 

-V 

1  •  u  «  M»  P  •  1  •  x  <  MP  =  J.  Therefore,  <T,u,J>  ,  as  defined  above, 
is  a  FSL  -  TOMA  scheme  under  which  f  is  feasible  and  J  <  P  •  (L-l) !  . 


Fact  4.8.  if  f  is  feasible  under  a  FSL-TDMA  scheme  <T,u,J>  ,  then 

.+ 

there  exists  a  TDMA  scheme  <T,x>  under  which  f  is  feasible. 

Proof.  By  letting  x  u  ,  we  obtain  the  desired  TDMA  scheme. 

"  J 

Thus,  we  have  shown  that  if  f  is  a  rational  vector,  then  feasibility 
under  TDMA  is  equivalent  to  feasibility  under  FSL-TDMA. 

Sometimes,  we  may  wish  to  find  a  FSL-TDMA  scheme  with  a  fixed  J. 

In  this  case  the  feasibility  of  f  is  equivalent  to  the  non-enqptiness 
of  the  region  defined  by  the  following  inequalities 

k  •  1 , . . . ,  L 

j  e  c.  i  i-i,...,j  i 

j  ■  1, . . . ,L  ;  k  ■  1, . * . ,L 

L 


•J 

E 

i-i 


- J  £i 


xij  +  xik  -  1 


xij  6  t0'1) 


i  m  1 , . . . ,  J  >  j  ■  I,.*., 


We  interpret  the  hihary  variables  * 


as  follows » 


1  if  a  transmission  over  link,  j 
is  to  be  made  at  slot  i 

0  otherwise 


If  the  above  region  is  empty  for  a  certain  J,  t  may  still  be 
feasible  under  FSL-TDMA  for  some  J'V  J* 

We  finish  this  chapter  by  stating  the  NP-completeness  of  a  restricted 

’ 

version  of  the  f-feasibility  problem  under  FSL-TDMA. 


Instance:  PRN  G  =  (N,A) ,  mean  link  arrival  rate  vector  f  such  that 

P  •  f  has  integer  elements  for  some  fixed  -  but  arbitrary  -  positive 

integer  P,  integer  J  such  that  3  <  J  <  P  •  (L-l) * 

o  —  o— 

Question:  Does  there  exist  a  FSL-TDMA  scheme  <T,u,J>  such  that  J  <  J 
and  T  u  >  J  f  ? 


Comment:  This  problem  remains  NP- complete  for  all  fixed  Jq  3.  (As 
stated  above,  JQ  is  part  of  the  input  and  is  arbitrary  in  the  specified 
range.) 

We  omit  the  proof  of  this  fact,  but  let  us  note  that  a  simple  proof 
is  possible  by  a  transformation  from  the  PARTITION  INTO  CLIQUES  problem 
(see  page  193  of  (GJ]  for  the  definition) . 


CHAPTER  V.  POLYNOMIAL  TIME  APPROXIMATION  ALGORITHMS  FOR  FF  AND  RF 


In  this  chapter  we  discuss  polynomial  time  approximation  algorithms 
for  the  feasibility  problems,  and  show  why  it  is  difficult  to  solve  them 
efficiently,  eyen  in  an  approximate  sense. 

5.1  Terminology  [PS] 

Let  P  be  an  optimization  problem  with  a  positive  cost  function; 

A 

let  I  be  an  instance  of  P;  let  V(I)  be  an  optimal  solution  to  I;  and 
let  c(V(I))  be  the  cost  of  0(1). 

Let  A  be  an  algorithm  which  returns  a  feasible  solution  V  (I) 

cl 

when  supplied  with  an  instance  1;  and  let  c(V  (I))  be  the  cost  of 

cl 

v  (I) . 

a 

The  algorithm  A  is  called  an  e  -  approximate  algorithm  if, 
for  any  instance  I  of  P,  we  have 

A 

|c(v  (I))  -  c(v(D)  1 
- S__ - <  e 

c(V(I) ) 

An  8  -  AA  (e-approximate  algorithm)  which  operates  in  polynomial 
time  is  called  an  e  -  PTAA  (6  -  polynomial  time  AA) . 

The  solution  V  (I)  returned  by  an  e-AA  A  is  called  an  6  -  approximate 

a 

solution. 

A  polynomial  time  approximation  scheme  (PTAS)  for  P  is  an 
algorithm  which,  when  supplied  with  an  instance  I  of  P  and  a  number 
6  >  0,  returns  an  6  -  approximate  solution  in  time  bounded  polynomially  . 
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in  the  size  of  the  instance,  where  the  polynomial  depends  on  6  . 

A  PTAS  is  called  a  FPTAS  (fully  PTAS)  if  its  running  time  is 
bounded  polynomial ly  in  1/e  and  the  size  of  the  instance. 

5.2  Negative  Results  About  FF 

The  optimization  version  of  FFt  which  we  shall  still  denote  by  FF, 

is  the  problem  of  finding  a  vector  x  >  ($,  such  that,  for  a  given  mean 

♦  +  "*  ■*  -► 

link  arrival  rate  vector  f,  Tx  >  f  and  1  •  x  is  minimized.  This 
problem  has  a  solution  whether  or  not  f  is  feasible.  We  shall  denote 
the  cost  of  an  optimal  solution  to  an  instance  I  »  <G,f>  of  FF  by 

Z (1) t  the  cost  of  an  approximate  solution  will  be  denoted  by  Z  (I). 

& 

We  now  list  the  main  results  about  FF  and  prove  them  in  the 
following  sections. 

Theorem  5.1. 

For  FF,  there  are  two  possibilities:  either  there  is  no  6  -  PTAA 
for  any  value  of  e,  or  there  exists  a  PTAS  for  all  8  >  0. 

Theorem  5.2. 

If  there  is  a  PTAS  for  FF,  then  there  is  a  PTAS  for  MC  (i.e.  the 
optimization  version  of  the  maximum  clique  problem) . 

Theorem  5.3. 

Unless  P  «  NP,  there  is  no  FPTAS  for  FF. 

We  are  unable  either  to  find  an  €  -  PTAA  for  FF  or  to  prove  that. 
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unless  P  “  NP,  there  exists  no  6  -  PTAA  for  FF.  However,  in  view  of 
Theorem  5.2  and  the  fact  that  a  PTAS  for  MC  -  an  extensively  studied 
problem  -  is  yet  to  be  found,  the  evidence  about  the  difficulty  of 
finding  an  6 -PTAA  for  FF  is  conclusive. 


5.3  Proof  of  Theorem  5.1 
•* 

Let  I  *»  <G,f>  be  an  instance  of  FF.  We  define  cxi  to  be  the 
instance  <G,  af>  where  a  >_  0. 

The  instance  I  **  1(1)  -  <G,f>  is  defined  as  follows: 


G  »  (N,A) 

N  **  {a^  :  a£N,  ie  {iGA  :  f±  >  o}} 

*«•  •*»>  ^ 

A  -  {(ai,bi)  :  a.jGN,  ^GN,  (afb)6A>U 

{(ai,bj),  (bj,ai)  :  i  t  j,  a.GN,  b^GN ,  iecj 


i  -  j 

i  i  j 


(Note  that  f^  is  the  mean  link  arrival  rate  on  link  i  of  G;  f^  is 
the  rate  on  link  (a,b)  of  G.) 


Since  the  case  f  ■  0  is  trivial,  we  shall  always  assume  that  there 
is  at  least  one  link  in  G  with  a  positive  mew  packet  arrival  rate, 
and  we  shall  let  k  denote  the  number  of  links  in  G  with  positive 
mean  arrival  rates.  The  links  will  be  labeled  such  that  f^>0, . . . ,f^>0. 


and  f. 


k+1 


0. 


4 
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The  subgraph 

Gl  -  ({a^SN  i  aSN),  i  (a,b)6A})  •  • 

—  th  % 

of  G  will  be  referred  to  as  the  i  copy  of  G  in  G,  because  G^  is  isomorphic 

to  G.  Thus,  there  are  k  copies  of  G  in  G. 

Before  we  proceed  with  the  proof,  let  us  illustrate  the  graph 

«■*  m*  •* 

G  “  (N,A)  by  an  example. 

Example.  Let  I  ■  <G,f>  be  as  in  Figure  5.1. a.  The  instance  I  »  <G,f> 
is  shown  in  Figure  S.l.b. 

The  cross-hatched  areas  between  the  copies  represent  links  ; 
thus,  e.g.,  each  node  in  copy  G^  is  connected  to  each  node  in  copy 
by  directed  links,  and  vice  versa. 

In  general,  for  any  two  copies  G^  and  G^ ,  there  are  two  cases: 

(1)  If  ieCj  in  G,  i.e.  if  link  iSA  conflicts  with  link  jGA,  then  each 

link  in  copy  G^  conflicts  with  each  link  in  G^j 

(2)  If  i0C^  in  G,  then  there  is  no  link  in  G^  that  conflicts  with  any 

link  in  Gj . 

The  mean  link  arrival  rate  vector  associated  with  the  links  in  G^ 

+  - 

will  be  denoted  by  f^.  The  instance  I  contains  k  instances  i^,...,!^ 
where  -  <G^,fj>  .  Notice  that  1^  and  f^I  are  identical. 

We  shall  first  argue  that  Z(I)  £  (Z<I))  .  Let 
V(I)  »  <u^,y^-j  i»i,...,|A|>  be  an  optimal  solution  to  I  where  u^ 
is  a  transmission  vector  of  G  and  y^  is  the  time  duration  associated 
with  u^. 
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The  system  £^V(I)  ■  <u^,  s  i  »  1, . . . , |a|>  is  an  optimal  solution 

A 

to  the  instance  f .1  and  has  a  cost  of  f^Z(Z)  for  each  j  **  l,...,k. 

A 

Thus,  V(I)  can  be  used  to  solve  each  of  the  instances  1^  (i=l,...,k) 

optimally.  Moreover,  the  solutions  of  these  instances  can  be  performed 

concurrently  since  certain  collections  of  copies  of  G  are  independent 

of  each  other.  To  be  more  precise,  let  u  be  a  transmission  vector  of  G 

-► 

such  that  u^  «  0  for  i  >  k,  i.e.,  u  uses  only  those  links  of  G  with  non¬ 
zero  mean  packet  arrival  rates.  Consider  the  set  {g^  :  ■  1}  ;  copies 

of  G  in  this  set  are  independent  in  the  sense  that,  if  s^  is  a  trans¬ 


mission  set  of  G. ,  then 

i 


IK 


is  a  transmission  set  of  G. 


ie{i:u^*l} 


It  is  now  easy  to  see  that  there  is  a  feasible  solution  to  I  with 
a  cost  of  (Z(I))2.  Therefore,  Z  (I)  <_  (Z(I))2. 

Now,  we  shall  give  a  reverse  argument.  A  feasible  solution  v (I) 
to  I  induces  a  solution  to  each  instance  1^  (i«l,...,k).  More 

precisely,  if  V(l)  -  <t^,x^  j  j  -  l,...,r>  ,  then  the  induced  solution 
V(I^)  to  1^  is  the  system  <u^,  y^  :  j  ■  l,...,r>  where  u^  is  the 
0-1  column  |a|  -  vector  obtained  by  deleting  all  rows  of  t^  that 
correspond  to  links  which  do  not  belong  to  and 


/  x. 


ri3 


-  < 


«u  ** 0 


otherwise. 


The  reader  should  verify  that  u. .  is  a  transmission  vector  of 

r 


G,  (i*l,...,k;  j»l,...,r),  and  j  u  .  y. .  >  f.  (i« 
1  j-1  13  1 


!1 , . . .  ,k)  • 
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r 

The  cost  of  V(I . ) ,  denoted  by  z. ,  is  ][  y.  .  . 

1  1  j=l  x:J 

Let  us  associate  a  column  0-1  |  A  |  -  vector  q_^  with  transmission 
■■ 

vector  t^  of  G  as  follows: 

1  je{l,...,k}  and  at  least  one  link  in  G. 

-*■  3 

is  used  by  t^. 

0  otherwise  . 

-► 

It  should  be  clear  is  a  transmission  vector  of  G. 

Mow,  consider  the  instance  I"  =  <G,f'>  where 


Z ^  i  * 

.0  i  =  k+1, ,  |a{ 


The  system  V(l')  =  <q^ ,x^  :  i  «  l,...,r>  is  a  feasible  solution 

to  l'  and  the  cost  of  V(l"),  denoted  by  Z(l'),  equals  Z(l). 

Let  us  define  a  =  min  {(Z^/f^):  i  =■  l,....,k}  ,  and  compare 

the  instance  I  \.ith  77  l'  .  Since  the  mean  link  arrival  rates  of 

a 

~  l'  dominate  the  corresponding  mean  link  arrival  rates  of  I,  the 

1  -*  Xi 

system  —  V(I  )  »  <q^,  —  :  i  ■  l,...,r>  is  a  feasible  solution  to  I 
with  a  cost  of  —  Z(l"). 

Therefore,  any  feasible  solution  V(i)  to  I  induces  k+1  feasible 

solutions  to  I;  k  of  which  are  on  the  copies  of  G,  and  one  is  on  the 

instance  —I*.  The  solution  induced  on  G,  is  the  system 
Q  1 

Vd^  ■  <*1^,  1  j“l, . . . ,r>  ,  and  has  a  cost  of  (Z^Vf^. 

The  induced  solution  with  the  least-cost  has  a  cost  of 
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Z (I)  ■  min  {a,  (Z(l')/o)}  since  a  «  min{ (Z, /f . )  :  i=l,...,k}.  Therefore, 
Z(z)  -  z(l")  >  Z(l)  a  >  (Z(I))2. 

The  optimal  costs  associated  with  Z  and  Z  must  also  satisfy 
E(i)  £  (Z(Z))2.  Thus  we  conclude  that  z(I)  «*  (z(I))2. 

Now  suppose  we  have  an  6  -  PTAA  for  FF;  when  we  apply  this 
algorithm  to  an  instance  Z  of  FF,  we  are  guaranteed  a  solution  V  (Z) 

zm  a 

such  that  - -  <  1  +  e  ,  where  Z  (Z)  =■  c(V  (Z)). 

Z(Z)  a  a 

Zf  we  apply  this  e  -  PTAA  to  Z  instead  of  Z,  we  still  obtain 

za(i) 

a  solution  st  h  that  x~z — £  1  +  e.  But  then,  by  searching  through 

Z(Z) 

the  induced  solutions,  we  can  obtain  a  solution  V(Z),  whose  cost 

satisfies  —  — ^  <  /l  +  6  .  Thus,  the  existence  of  an  e  -  PTAA 

Z(Z) 

implies  the  existence  of  a  (  /I  +  €  -  1)  -  PTAA.  (We  leave  it  to  the 

reader  to  prove  that  the  new  algorithm  is  still  polynomial  time.) 

Suppose,  for  an  instance  Z  of  FF,  a  5-  approximate  solution  is 

desired  (<5>0) .  How  do  we  use  the  6  -  PTAA  -  which  we  assume  exists  - 

to  obtain  the  desired  accvrracy?  Zf  6  e,  then  there  is  no  problem. 

Let  us  consider  the  case  6  <  e,  and  let  m  be  the  smallest  integer 

2~m 

such  that  1  +  {  >  (1  +  6)  .  Then,  we  apply  the  6  -  PTAA  to  the 

mth  order  instance  Z  ^ ,  which  is  obtained  from  Z  by  recursively 

mm 

applying  the  construction  of  Z  m  times,  and  descend  to  a 
2“* 

1(1  +  8)  -  1)  -  approximate  solution  tc  Z  by  following  the 

least-cost  induced  solutions  at  each  stage. 

This  con^letes  the  proof  of  Theorem  5.1. 
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5.4  Proof  of  Theorem  5.2 

The  proof  is  based  on  two  observations : 


(1)  There  exists  a  PTAS  for  MTS  if  there  exists  a  PTAS  for  FF, 

(2)  There  exists  a  PTAS  for  MC  if  there  exists  a  PTAS  for  MTS. 


The  second  observation  readily  follows  from  the  transformation 
algorithm  Al  of  Chapter  IV;  we  leave  out  the  proof. 

TO  prove  the  first  observation,  we  shall  consider  the  transform¬ 
ation  algorithm  A2  of  Chapter  IV. 

Suppose  there  is  a  PTAS  for  FF  and  we  wish  to  obtain  a 
y  -  approximate  solution  to  an  instance  I  «*  <G,B>  of  MTS,  where 
G  »  (N,A)  is  a  PRN,  and  B  £L  A.  By  using  A2,  we  can  transform  I  to 
an  instance  of  FF.  But,  there  is  a  difficulty  since  we  do  not  know 
what  value  of  k  to  use  in  the  transf ormation .  The  difficulty  is 
overcome  by  trying  all  values  of  k  between  1  and  | B | .  For  each 
value  of  k,  we  transform  the  instance  I  of  NTS  to  an  instance  of  FF, 
and  we  operate  the  PTAS  for  FF  with  an  accuracy  requirement  of  6  such 


that  0  <  6  <  '(i  y)  (l"  |b|  +  |b|2)  *  Eac^  tiro®  we  obtain  an  approximate 

solution  to  MTS  by  simply  searching  through  the  transmission  vectors 

used  by  the  PTAS.  We  claim  that  the  best  of  these  approximate 

solutions  is  a  y-  approximate  solution  to  the  instance  I  of  MTS. 

* 

To  see  why,  suppose  the  optimal  value  of  the  instance  I  is  k  .  When 

* 

the  PTAS  is  operated  on  the  instance  of  FF  with  k  «•  k  ,  an 


e  -  approximate  solution  must  have  a  cost  less  than  or  equal  to 
1  +  6,  because  the  optimal  cost  equals  1.  But  any  feasible  solution 
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must  spend  at  least  6k  ( ] B [  -  1)  fraction  of  time  in  satisfying  the 
service  requirements  of  the  links  in  D.  In  the  remaining 

(1  +  e)  -  6k  ( I B I  -  1)  =6+6  fraction  of  time  the  PTAS  must  use 

* 

transmission  vectors  which  on  the  average  enable  at  least  6k  / (6+6) 
v 

links  of  B  .  Therefore,  the  PTAS  uses  a  transmission  set  which 

constitutes  a  e/(6+6)  -  approximate  solution  to  the  instance  I  of 

6  Y 

MTS.  If  we  solve  £  y  for  6,  we  obtain  B  £  — y  ^  •  Replacing  6 

by  its  smallest  value  as  k  ranges  from  i  to  |b|,  i.e.  by  '£')  • 

we  see  that 


e  < 


-  (l-Y)  (1-  |  B |  +  jB  j 


guarantees  ay-  approximate  solution  to  the  instance  I  of  MTS. 


5.5  Proof  of  Theorem  5.3 


In  section  5.4,  we  have  seen  that  a  PTAS  for  FF  can  be  used  to 
obtain  approximate  solutions  to  instances  of  MTS. 

In  this  section  we  shall  prove  that  if  there  is  a  FPTAS  for  FF, 
then  we  can  obtain  an  optimal  solution  to  any  instance  of  MTS  in 
polynomial  time. 

* 


Suppose  1=  <G,B>  is  an  instance  of  MTS.  Let  k  denote  t>.e 
optimal  value  pf  I. 

As  in  the  previous  section,  we  shall  transform  I  to  an  instance 
of  FF  for  all  values  of  k  between  1  and  [b|.  When  k»k  ,  the  FPTAS 
for  FF,  which  we  assume  exists,  guarantees  a  solution  V(I)  such 


that 


c(V(I) )  > 
k* 


6 _ 

6  +  e 


where 


6 


We  would  like  to  choose  6  small  enough  so  that  - ? —  <  -  , 

k  6  +  6 

* 

because  this  would  guarantee  that  c(v(X))  ■  k  . 

Solving  for  e,  we  obtain  the  following  condition: 

|  >(k*  -  1) (1  +  k*( |b|  -  1)). 

* 

Since  we  do  not  know  the  value  of  k  in  advance,  we  choose 

|  -  (|b|  -  1) <1  +  Jb|  (|b|  -  1))  +  1  . 

We  note  that  ^  required  for  obtaining  am  optimal  solution  is 
polynomially  bounded  in  the  size  of  the  problem.  This  completes  the 
proof. 

5.6  Negative  Results  About  RF 
Theorem  5,4 

There  exists  a  PTAS  for  FF  if  there  exists  a  PTAS  for  RF. 

There  exists  an  6  -  PTAA  for  FF  if  there  exists  an  6-PTAA  for  RF. 

Unless  P  «  NP,  there  exists  no  FPTAS  for  RF. 

The  proof  of  this  theorem  follows  immediately  when  one  considers 
the  transformation  algorithm  A3  of  Chapter  IV. 

The  implications  of  Theorem  5.4  about  the  existence  of  approximation 
algorithms for  the  r- feasibility  problem  should  be  cleau:  in  view  of 


Theorem  5.2 
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CHAPTER  VI.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 

At  the  end  of  this  thesis,  we  can  say  with  some  certainty  that 
PRNs  are  complex?  even  the  most  fundamental  problems  are  not  likely 
to  be  solved  in  a  practical  way.  For  example,  if  NP  /  P,  there 
exists  no  polynomial  time  algorithm  for  deciding  whether  a  given  set 
of  mean  o-d  arrival  rates  belongs  to  the  capacity  region  of  a  given  PRN. 

We  have  proved  several  negative  results  about  TDMA  schemes;  let 
us  note  as  an  extension  that  similar  results  are  readily  obtained 
for  FDMA  (frequency-division-multi-access)  schemes  by  using  a  PRN 
model  in  which  each  station  is  capable  of  transmitting  and  receiving 
simultaneously  as  many  packets  as  desired,  provided,  of  course, 
the  frequency  bands  that  are  used  do  not  overlap.  We  should  remind 
the  reader  that  the  above  conclusions  are  based  on  worst-case  considera¬ 
tions  and  should  be  interpreted  accordingly. 

Our  initial  interest  in  studying  TDMA  schemes  in  the  course  of 
this  thesis  work  was  to  obtain  a  benchmark  against  which  the  perform¬ 
ance  of  other  multi-access  schemes  could  be  compared.  As  is  implicit 
in  the  proof  of  Theorem  3.2,  no  multi-access  scheme  is  superior  to 
TDMA  as  far  as  throughput  is  concerned.  However,  there  are  situations 
where  TDMA  is  undesirable,  because  it  leads  to  extremely  long  delays? 
eg.  a  single-receiver  PRN  with  many  small  users  is  one  such  case. 

When  TDMA  is  undesirable,  a  multi-access  scheme  such  as  Aloha  becomes 
attractive.  In  this  sense  TDMA  and  Aloha  are  complementary  to  each 
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other,  but  this  tradeoff  still  remains  to  be  quantified. 

As  a  second  suggestion  for  research,  TDMA  schemes  . can  be  studied 
in  PRNs  where  connectivity  between  the  nodes  of  the  network  depends 
on  the  distance  between  them. 


y 
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